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O ■ 1 Introduction 
(N , 

The action of genes is manifested statistically in sufficiently large communities of 
matching individuals (belonging to the same species). These communities are called 
^ I populations [2]. The population exists not only in space but also in time, i.e. it has 
H ' its own life cycle. The basis for this phenomenon is reproduction by mating. Mating 
- - - in a population can be free or subject to certain restrictions. 

The whole population in space and time comprises discrete generations -Fq) -^ij 
The generation Fn+i is the set of individuals whose parents belong to the gen- 
eration. A state of a population is a distribution of probabilities of the different 
types of organisms in every generation. Type partition is called differentiation. The 
simplest example is sex differentiation. In bisexual population any kind of differ- 
entiation must agree with the sex differentiation, i.e. all the organisms of one type 
must belong to the same sex. Thus, it is possible to speak of male and female types. 

The evolution (or dynamics) of a population comprises a determined change 
of state in the next generations as a result of reproductions and selection. This 
evolution of a population can be studied by a dynamical system (iterations) of a 
quadratic stochastic operator. 



1 



The history of the quadratic stochastic operators can be traced back to the 
work of S. Bernshtein [1]. For more than 80 years this theory has been developed 
and many papers were pubHshed (see fT]-[15]). Several problems of physical 
and biological systems lead to necessity of study the asymptotic behavior of the 
trajectories of quadratic stochastic operators. 

Let E = {1,2, ...,m}. By the (m — 1)— simplex we mean the set 

m 

= {x = ixi,...,xm) e i?"^ : X, > 0,^x, = 1}. (1) 

1=1 

Each element x E S"""^ is a probability measure on E and so it may be looked 
upon as the state of a biological (physical and so on) system of m elements. 

A quadratic stochastic operator V : 5""*"^ — > 5"""^ has the form 

m 

^ ■ ^'k = ^ Pij,kXiXj, (fc = l,...,m), (2) 
where Pij,k— coefficient of heredity and 

m 

Pij,k = Pji,k>0, '^Pij,k = '^, iij,k = l,...,m). (3) 

A;=l 

For a given x^^^ € S"^~^, the trajectory {x^"'^}, n = 0, 1,2, ... of x^'^^ under the 
action of QSO ([2]) is defined by x("+i) = y(a;(")), where n = 0, 1, 2, ... 

One of the main problems in mathematical biology is to study the asymptotic 
behavior of the trajectories. There are many papers devoted to study of the 
evolution of the free population, i.e. to study of dynamical system generated 
by quadratic stochastic operator ([2]), see e.g. [3]-[ll]. In [13] a survey of theory 
quadratic stochastic operators is given. 

In this paper we find a condition under which the evolutionary operators of 
bisexual population is contraction. 



2 Definitions 

In this section following [2], we describe the evolution operator of a bisexual popu- 
lation. Assuming that the population is bisexual we suppose that the set of females 
can be partitioned into finitely many different types indexed by {1, 2, n} and, sim- 
ilarly, that the male types are indexed by {1, 2, v}. The number n -\- v is called 
the dimension of the population. The population is described by its state vector 
(x, y) in 5""^^ x S'^~^, the product of two unit simplexes in and M*^ respectively. 
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Vectors x and y are the probability distributions of the females and males over the 
possible types: 

n u 

> 0, ^Xi = 1; Vj > 0, ^ Vj = 1. (4) 

i=l 3=1 

Denote S = S"'~^ x S'^~^. We call the partition into types hereditary if for 
each possible state z = (x,y) e 5 describing the current generation, the state 
z' = {x', y') G 5 is uniquely defined describing the next generation. This means that 
the association z ^ z' defined a map V : S ^ S called the evolution operator. For 
any point z^^^ € S the sequence z^^^ = V{z^^^^^),t = 1,2, ... is called the trajectory 
of z^^"^ . Let p[j \ and p^J^^ be inheritance coefficients defined as the probability that 
a female offspring is type k and, respectively, that a male offspring is type /, when 
the parental pair is k = 1,2, n; and j, I = 1,2, u). We have 



k=l 1=1 



Let z' = {x', y') be the state of the offspring population at the birth stage. This 
is obtained from inheritance coefficients as 



W : < 



x'k = E Plik^iVj^ {l<k<n) 



(/) 

^iVj, \± :^ 1^ :^ It) 

(6) 



y'l = E P'i^ixiyj, {l<l<v) 
«J=1 



We see from ([6]) that for a bisexual population the evolution operator is a 
quadratic mapping of 5 into itself. But for free population the operator is quadratic 
mapping of the simplex into itself given by ([2]). 

In [8j an algebra of the bisexual population is defined as the following: 
Consider {ei, e„+i/} the canonical basis on M"+'' and divide the basis as = 

ei,i = l,...,n and e\"^^ = en+i,i = l,...,y. Introduce on M""'"'' a multiplication 
defined by 



M)Am) ^ (m) (f) _lf ^ Jf) Jf) ^ ^ Jm)Jm) 



if ^ if) if) , ^ M 

\ k=l 1=1 J 

ePei^)=0, i,k = l,...,n; (7) 

(m) (m) „ . , -I 

e} 'e\ - 0, 2,1 = l,...,v; 

Thus the coefficients of bisexual inheritance is the structure constants of an 
algebra, i.e. a bilinear mapping of M"^'^ x W^^" to M""*"*^. The general formula for 
the multiplication is the extension of d?]) by bilinearity, i.e. for z, t G 



n V n V 

i=l j=l i=l j=l 
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using ([7]), we obtain 

k=\ \%=\j=\ J /„N 

1=1 \i=lj=l J 

From ([8|) and using ([6]), in the particular case that z = t, i.e. x = u and y = v, 
we obtain 

= = E ( E E^^S/i^^^/,- ) ei^h 
k=i ^ i=i j=i ' 

+E(EE^'SrM^.)-^^ = ^(-)- 

for any z € 5. This algebraic interpretation is very useful. For example, a bisexual 
population state z = {x, y) is an equilibrium (fixed point) precisely when z is an 
idempotent element of the set S, i.e. z = z'^. 

The algebra B = By/ generated by the evolution operator W (see ([6|)) is called 
the evolution algebra of the bisexual population. 

In [8] it was shown that if z is a fixed point then z € R^'^'^ |J RJ^'^'^ , where 

n V 

i?:;+'^ = {z = (x,y): j;x, = J^y,=7?}, r/ = 0, 1. (9) 
For simplex 5 = 5"^^ x S^~^ by tangent space we get 

n V 

R^^-" = {z = {x, y) : E = E = (10) 
i=i j=i 



3 Contracting operators 

In operator theory, a bounded operator W : X ^ Y between normed vector spaces 
X and Y is said to be a contraction if its operator norm ||W^|| < 1. 

An extremal example of a quadratic contraction is the constant operator. In this 

( f) (m) 

case the coefficients 1, p,-^- / do not depend on i and j. This suggests that for a 
sufficiently small scattering of coefficient for every fixed k, I the quadratic operator 
will be a contraction. This remark can be expressed as a precise theorem. 
The Lipschitz constant of an operator W : M""'"'^ — >■ M"'+'' is 

\\Wz-Wt\\ 

L{W) = sup — . — , 

z^t \\z ~ i^W 

where || • || is some norm in M""'"'^. If this norm can be chosen so that L{W) < 1 then 
W will be a strict contraction in this norm with the consequences: unique fixed point, 
convergence of all trajectories to this point, exponential rate of convergence. Unless 
otherwise specified, we will use the h— norm in the basis = ei,i = l,...,n 
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and = en+i,i = defined as = + X] ^ — i^^v) 

i=i j=i 

n 



4 = 1 



J = l 



Lemma 1. f^. Let A be a convex n— dimensional compact in M", F : A — t- A be 
a smooth map. Then (for any norm) L(F) = max 

Lemma 2. f^. Let a matrix A = {aij)fj^i satisfies 

n n n 

^a,i = J^a,2 = ... 

Then 



0>in. . 



(11) 



1=1 i=l 



1=1 



1 

ll^ll^o II = 9 ™ax ^ laj,;, - a 



1=1 



(12) 



where A\&q is restriction operator A on Mq. 

For each z £ B we have hnear operator : B B defined by Mz{t) = zt. 
Theorem 1. The following inequality holds for the Lipschitz's constant 

L(W) < max ( Y. Iflfi. " f 'S,! + E l^.'S " ^Sil + 

''''^^fc=l Z=l ^ 

( X/ l^'iii,fc ~ -Piji.fcl + X/ l^ii!« ~ ^i2!il ) ■ 
^fc=l «=1 ^ 



max 



(13) 



Proof. For the operator in S the derivative is 



^ (/) 

Elm) 
PijM 

i=i 



\ i=i 



E PnjliVj E p. 



jl,l-^« 



E P£l^^ ^ 



j=l 



i=l 



E Pni,nyj E Pil,n^i 



E 

i=l 
n , , 
v-^ (m) 

E pI,i^^ 

i=l 



V 



n 



E(m) ^ (m) 

Ki/iyj E Pii> 

i=i 



Pnj,>j E P 



(m) 
il,i/^« 



Elm) 
pI,uX 

i=l 



(14) 



Am) 



k=l 



1=1 



where Mj/^ = M (/) and m/"^'* = M (m) is the multiphcation maps with matrixes 
(:p!il)^fc=i ^"^^ respectively (pj™^)^;^^. 
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By Lemmanwe have L{W) = 2max ||M^|| < 2max \\Mjf'\\ + 2max ||M^^' 
By Lemma m 

1 / " 



A:=l «=1 



□ 



Corollary 1. An evolutionary operator ^ is a strict contraction if 



--(EbS.-^S5.i+Ei^!S-^Si)+ 
jj- (EbS£.-^S£.i+Ei^SS-^Si) <i (15) 

For evolutionary operators with positive coefficients there is a multiplicative 
estimate of the distance from the evolutionary operator to the constant one. Let 

„{/) „M 
= fi^iW) = max ^™ = = max 

n,i2,j,k Ml 1,31,32,1 nj^"^', 

^i23,k ^^32,1 

and let C{W) equal to LHS of (15). 
Lemma 3. 

< 4 V4 + 7- (16) 

Proof. If a, P > and ^ = max(^, ^) then obviously 

|a-/3| = ^(a + /3). 
Ai + 1 

Hence 

I (/) _ (/) I ^ 1!L_ZJl( (/) _L ^ I (™) _ M I ^ ^"^ - 1 ^ (m) (m)x 

and respectively 

I (/) _ (/) I ^ /^-^ - 1 ^ (/) , (/) \ I M _ (m) I . ^"^ - 1 ^ (m) (m)N 
f'ija.fcl - ^/ _^ X '^Pijuk '^Pij2,k)' \Piji,l Pi32,l\ - i ^Piji,l ^ Pij2,l) ' 

It remains to sum these inequalities over k and respectively over I, keeping in 
mind that 

n n u 1/ 

Z^Piij,k — Z^Pi23,k ~ Z^Pijul ~ Z^Pij2,l ~ 
k=l k=l 1=1 1=1 

□ 
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Corollary 2. 



Corollary 3. // Tfi^fi 
strict contraction. 

Corollary 4. Let fi 



L{W) < 4 



/ 



/ + 1 



1 



1 



/i"^ + 1 



(/i-^ + /U™') < 9 i/ien the evolutionary operator ^ is a 

max 



L(W^) < 8 



/x + 1 

anrf if fi < J then the evolutionary operator ^ is a strict contraction. 

Let us give several examples and check the condition of Corollary [TJ 
Example 1. Consider the operator 



W : < 



X 



X2 = jxiyi + ^xiy2 + lx2yi + 7X2^2, 
fXiyi + ^xiy2 + ^X2yi + 7X2^2, 
fXiyi + |xiy2 + ^^22/1 + f X2y2- 



2/1 



(17) 



The coefficients of the operator (jl7p as the following 



Pll,l — 


3 
7 


^'12,1 — 


1 
2 


^(/) _ 
P2l,l — 


1 
2 


Jf) 
P22,l 


M) _ 
^•11, 2 — 


4 
7 


„(/) _ 


1 
2 


Jf) _ 
/^21,2 — 


1 
2 


P22,2 


(m) 


4 
7 


(m) 
Pl2,l = 


1 
2 


(m) 
P2l'l = 


1 

2 


(m) 
^^22,1 


(m) 


3 
7 


(»n) 
^'12, 2 ~ 


1 
2 


(m) 
^'21, 2 ~ 


1 
2 


(m) 
^22,2 



It is easy to check that condition (fT5|) satisfied for (fT7|) . Indeed 

• n 
■ A:=l 



max 



(/) 

iij,k 



.Jf) I 



(m) 



1 , \ " I (m) (m) A 
;— 1 / 



Consequently, this operator is a strict contraction and it has unique fixed point 
^, |, ^). Moreover any trajectory of (fT7|) converges to the fixed point. 
The following example shows that the condition of Corollary [1] is not satisfied 
and evolutionary operator has periodic trajectory. 
Example 2. Consider the operator 



W : < 



y'l 
y'2 



■■ xiyi 

Xiy2 + X2 
X2y2 

xi + X2yi 



(18) 
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It easy to check that operator does not satisfy the condition of Corollary [TJ 
We rewrite the operator (jlSp in the form 



W : < 



Denote x„ = x^^\ yn = y^^ then from (fT9|) we have 

yn+1 = (1 - Xn)il - yn) 

Since < x„y„ < Xn from the first equation of ([2U|) it follows that lim = 

n— >oo 

X* = 0. Indeed, for (x°,y°) G mt(5^ x S^) we get from (|20]l 

^ = a-xjfl-^ 



yi = (l-xi)(l-yi) 



Xji-^-l \ Xf) 

3^n+23^n — ^n+l(l 2;j2)(Xn ^Jj^-j-i), 

lim x„+2X„ = lim x„+i(l - x„)(x„ - x„+i), 

ra— ^oo n— ^-oo 

{x*f = 0, X* = 0. 

Now consider the operator 

x' = xy — x'^y — xy^ + x^y^ 
: I , 2 2 2 2 (21) 

y = X + y — xy — x y — xy + x 7/ 

Clearly, the operator W'^ has fixed points (0, y), < y < 1. The point (0, y) is a 
saddle point. 

It is easy to check that the set {(x,y) ^ x : xi = 0} is an invariant 
subset for (jlSp . Any point of the invariant subset is periodic point with period 
two for operator (llSp . So trajectory of the operator with an initial point form 
invariant subset does not converge. Thus operator (jlSp has a trajectory which does 
non-converge to the fixed point (0, 1, ^, |). 

The following example shows that condition of Corollary [1] is sufficient but is not 
necessary. 

Example 3. Consider the operator with coefficients of inheritance 



„{/) 


= 


„(/) _ 

^•12,1 — 





„{/) 
P21,l 


1 

~ 2 


„(/) 

P22,l 


Jf) 

Pu,2 


= 1 


Jf) _ 

Pl2,2 — 


1 


Jf) 
P21,2 


1 
2 


P'22,2 


(m) 

Pui 


= 


(m) 


1 

2 


(m) 
^21,1 


= 


(m) 
^^22,1 


(in) 


= 1 


(rra) 
^•12,2 — 


1 
2 


(m) 
^21,2 


= 1 


(m) 
P22,2 
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i.e. the evolution operator has the form 



W : < 



x'2= Xi + \X2 

, 1 (22) 

Vl = 22/2 



2/2=2/1 + \V2 

It easy to check that operator does not satisfy the condition of Corollary [TJ 



(n V 
fc=i 1=1 

, n V ^ 

( E - fS£,J + E - f.'Si) = 2 > 1. 

But any trajectory of ([^^ converges to (i, |, i, |). 



Indeed, from ()22p we have 

.r^)=i(i-x(")) 

We consider following one dimensional dynamical system. 

/(x) = ^(l-x) 

It has unique fixed point x = \^ and decreasing on [0,1]. Easy to check that 
j'{x) = —\ and |/'(|)| = \ < 1 therefore the fixed point x = ^ is attracting. 

We claim that any trajectory of f{x) converges to the fixed point x = ^. Indeed, 
we have 

k=l 

and 

1- ,2n, , r A 22" - 1 x \ 1 
hm f (x) = lim - • K h ^ = -, 

lim = lim ( - ■ + ^i— - = -. 

n-s-oo'' ^ ^ n^oo\^3 22" ^ 22"+l 22"+ly 3 

So for any initial point trajectory of (f22]) converges to (|, |, ^, |). 
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